Chiral charge order in lT-TiSe2: importance of lattice degrees of freedom 



m 



o 


12 


(N 


13 




14 


u 




Ph 




< 


17 


l> 


IS 




19 




20 


<D 


21 


1 


22 


;h 





t« 



c^ 



X3 
O 

o 



> 

in 

o 
m 



X 



B. Zcnkcr\ H. Fehske\ H. Beck^, C. Monney^ and A. R. Bishop'' 

^InsUtut fiir Physik, Ernst-Moritz-Arndt-Universitdt Greifswald, D-17489 Greifswald, Germany 

^ Departement de Physique and Frihourg Center for Nanomaterials, 

Universite de Frihourg, CH-1700 Fribourg, Switzerland 

'^ Fritz- Haber-Institut der Max Planck Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany 

* Theory, Simulation, and Computation Directorate, 

Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA 

(Dated: April 9, 2013) 

We address the question of the origin of the recently discovered chiral property of the charge- 
density-wave phase in lT-TiSe2 which so far lacks a microscopic understanding. We argue that the 
lattice degrees of freedom seems to be crucial for this novel phenomenon. We motivate a theoretical 
model that takes into account one valence and three conduction bands, a strongly screened Coulomb 
interaction between the electrons, as well as the coupling of the electrons to a transverse optical 
phonon mode. The Falicov-Kimball model extended in this way possesses a charge-density-wave 
state at low temperatures, which is accompanied by a periodic lattice distortion. The charge ordering 
is driven by a lattice deformation and electron-hole pairing (excitonic) instability in combination. 
We show that in addition an explicit phonon-phonon interaction must be taken into account to 
achieve a stable chiral charge order. The chiral property is exhibited in the ionic displacements. 
Exploring the interplay between electron-electron interaction and the electron-phonon coupling in 
the whole parameter regime, we provide the ground-state phase diagram of the model and give an 
estimate of the electron-electron and electron-phonon interaction constants for lT-TiSe2 . 

PACS numbers: 71.45.Lr, 71.27.-|-a, 63.20. kk, 71.35.Lk, 71.38.-k 



I. MOTIVATION 

Charge-density-waves (CDWs) brought about by 
electron-phonon^ or electron-electron^ interactions are 
broken-symmetry ground states, typically of low- 
dimensional (D) solids with anisotropic properties.'^ A 
prominent material exhibiting such a periodic real-space 
modulation of its charge density is the transition-metal 
dichalcogcnidc ir-TiSc2. This quasi-2D system under- 
goes a structural phase transition at about 200 K, at 
which a commensurate 2x2x2 superstructure accom- 
panied by a CDW develops.'' Thereby the CDW features 
three coexisting components and, for this reason, is de- 
noted as triple CDW. Although the charge ordered phase 
in ir-TiSc2 has been a matter of intensive research for 
more than three decades, the driving force behind the 
phase transition has not been identified conclusively. 

Recent experiments on lT-TiSe2, pointing to a very 
unusual chiral property of the CDW, have reinforced the 
interest in this problem. ^'^ An object exhibits chirality 
if it cannot be mapped on its mirror image solely by 
rotations and translations. For a CDW phase character- 
ized by a scalar quantity such chirality has not been de- 
tected before. The scanning tunneling microscopy mea- 
surements performed by Ishioka and co-workers, how- 
ever, show that the amplitude of the tunneling current 
modulates differently along the CDW unit vectors in IT- 
TiSe2.^ Since the tunneling-current amplitude directly 
measures the local electron density, the charge density 
modulates differently along the three unit vectors. As 
a result the material in its low-temperature phase will 
not exhibit a three-fold symmetry as suggested by the 
triangular lattice structure. The Fourier transformation 



56 of the scanning tunneling microscopy data demonstrates 

57 a triple CDW as well as a different charge modulation 

58 along each CDW component with the respective order- 

59 ing vector Qa, a = 1, 2, 3.^ If one orders the triple-CDW 

60 components according to their charge modulation ampli- 

61 tude in ascending order, in a sense a direction is induced 

62 and the triple CDW exhibits chirality because the mirror 

63 symmetry is broken, in contrast to usual CDWs;^ see the 

65 schematic representation by Fig. 1. Note that clockwise 

66 and anticlockwise chiral CDWs were found in the same 

67 sample, suggesting that these states are degenerate. This 




FIG. 1. (color online) Electron-density pattern for a trian- 
gular lattice in case of (a) a nonchiral CDW or (b,c) chiral 
CDWs. Filled circles picture the value of the charge densities, 
where equal colors mark equal densities. For the nonchiral 
CDW shown in (a) the density modulation along the order- 
ing vectors Qi, Q2, and Q3 is equal. Reflection along an 
ordering vector yields the same density pattern, i.e., mirror 
symmetry exists. The situation changes for a chiral CDW. 
A clockwise CDW (red arrow) is illustrated in (b). Now re- 
flexion along Qi yields the situation depicted in panel (c). 
Obviously the pattern (c) describes an anticlockwise CDW: 
That is, for a chiral CDW mirror symmetry is broken. 



68 two-fold symmetry is corroborated by optical polarime- 

69 try measurenients.^ Ishioka and co-workers furtherraore 

70 noticed that the experimental data can be reproduced by 

71 a charge density modulation of the form 



ni{Qa) = ^COs(QaRi + 6a), 



(1) 



where A is the modulation amplitude and 9a are initial 
phases.^ For a chiral CDW to exist the phases of the 
CDW components must differ, i.e., Qx^Q-i^^^z- 

From a theoretical point of view the chiral CDW in IT- 
TiSe2 was addressed by a Landau- Ginzburg approach. ^'^ 
Thereby the relative phases of the CDW order parame- 
ters were obtained by minimizing the free energy func- 
tional. Two CDW transitions were found with decreas- 
ing temperature: Firstly a standard (nonchiral) CDW 
appears, and subsequently a chiral CDW emerges, i.e., 
TnonchiraiCDW > TlhiraiCDW- Within the CDW phasc 
three distinct orbital sectors are occupied, leading to an 
orbital-ordered state and three interacting lattice dis- 
placement waves (with different polarizations). 

An open issue is the microscopic mechanism driv- 
ing the CDW transition. Basically two scenarios were 
discussed in the literature, where the charge order re- 
sults from purely electronic, respectively electron-lattice, 
correlations. Angle-resolved photoemission spectroscopy 
data reveal a relatively large transfer of spectral weight 
from the original bands to the back-folded bands (due to 
the CDW transition) , compared with the small ionic dis- 
placement. This suggests an electronic mechanism within 
the excitonic insulator (EI) scenario. ^'^° A corresponding 
tight-binding calculation estimates the amplitude of the 
lattice deformation caused by an EI instability to be of 
the same order as the measured one."'^"'^ The gradual sup- 
pression of the CDW phase by changing solely electronic 
properties by intercalation with S or Te further corrobo- 
rates the EI concept. ^^ Most convincingly, time-resolved 
photoemission spectroscopy reveals an extremely fast re- 
sponse of the CDW to external light pulses which favors 
an electronic mechanism. ^'^ Alternatively, the coupling 
to the lattice degrees of freedom may drive the CDW 
transition, e.g., by a cooperative Jahn- Teller effect. ^^'^^ 
Here the particular form of the phonon dispersion and the 
softening of transverse optical phonon modes were elabo- 
rated within a tight-binding approach and found to agree 
with the experimental results. ^^"^^ The same holds for 
an ab-initio approach^° to a Jahn- Teller effect. Likewise 
the onset of superconductivity by applying pressure may 
be understood within a phonon-driven CDW scenario. ^^ 
Since some properties of the CDW in ir-TiSe2 can be 
understood by the excitonic condensation of electron- hole 
pairs and others by the instability of a phonon mode, a 
combined scenario has been proposed. ^^ 

As yet it is unclear whether the chiral property of the 
CDW favors the electronic or lattice scenario, or a com- 
bination of both. In the present work, this issue is ad- 
dressed amongst others. We start by investigating the 
CDW from an EI perspective. The corresponding mean- 
field approach for an extended Falicov-Kimball model is 



124 presented in Sec. II A. We show that the EI scenario is 

125 insufficient to explain a stable chiral CDW. We proceed 

126 by including the lattice degrees of freedom. It turns out 

127 that, besides an electron-phonon coupling, the plionon- 

128 phonon interaction has to be taken into account in order 

129 to stabilize chiral charge order. This is elaborated in 

130 Sec. II B. The self-consistency equations for the CDW 

131 and EI order parameters are derived in Sec. II C. The 

132 CDW state is characterized analytically in Sec. II D. Scc- 

133 tion III contains our main (numerical) results. Here we 

134 give the functional dependences of the EI order param- 

135 eter and the lattice distortion on the model parameters 

136 and temperature, derive the ground-state phase diagram, 

137 and estimate the interaction constants for lT-TiSe2. In 

138 Sec. IV. wc summarize and conclude. 



II. MODEL AND THEORETICAL APPROACH 

A. Electronic degrees of freedom 

1. Band structure 

Since the electronic properties of lT-TiSe2 are dom- 
inated by the electrons near the Fermi energy, in what 
follows we only take into account the top valence band 
and the lower conduction band. The maximum of the 
valence band is located at the F-point. The conduction 
band exhibits minima at the three L-points, see Fig. 2. 
To facilitate the notation, we artificially split the conduc- 





FIG. 2. (color online) First Brillouin zone (BZ) of ir-TiSe2 
with high symmetry points in the normal phase (solid line) 
and in the CDW phase (dot-dashed line). Red arrows show 
the CDW ordering vectors. Left panel: projection onto the 
xj/-plane, right panel: projection onto the yz-plane. 

150 

151 tion band into three symmetry-equivalent bands indexed 

152 by a, each having one minimum at the point La- The 

153 band dispersions of these three conduction bands mimics 

154 the true band structure close to the L-points.^° Figure 3 

156 illustrates the situation close to the Fermi level. Then 

157 the free electron part is written as 



Hc = ^ewfflfi, + ^ei,a< 



■^ c 



(2) 



k,a 




(a)band structure along (b)band dispersion close to 

high-symmetry directions of the BZ the Fermi level 

FIG. 3. (color online) Model band structure in the normal 
phase. The valence band is colored in black and the conduc- 
tion bands are colored red, blue, and green. 



158 where /]^ annihilates (creates) an electron in the valence 

159 band with momentum k and c^^ annihilates (creates) an 

160 electron in the conduction band with momentum k and 

161 band index a. The corresponding valence-band disper- 

162 sion and the conduction-band dispersions are denoted as 

163 eic/ and £ka- They will be specified in Sec. Ill A. The 

164 spin of the electrons is neglected. 

165 Taking the band structure and a band filling factor 

166 n = 1/4 into account, ir-TiSe2 resides in the vicinity of 

167 a semimetal-semiconductor transition, cf. Fig 3. Accord- 

168 ingly the chemical potential jjl is determined by 



V 



E^ 



1, 



wher 



k k 

*" = 1^ E«) = W E(4aCkc> 



(3) 

(4) 
(5) 



180 where Ucc denotes the Coulomb repulsion among the con- 

181 duction electrons. The on-site Coulomb interaction Ufc 

182 between valence and conduction band electrons deter- 

183 mines the distribution of electrons between these "sub- 

184 systems" and therefore may drive a valence transition, as 

185 observed, e.g., in heavy fermion and intermediate- valence 

186 Tm[Se,Te] compounds. ^"''^^ If the total electronic model 

187 contains an explicit hybridization between / and c clcc- 

188 trons^®'^^ or, as in our case, dispersive c and f bands, ^^ 

189 coherence between c and / particles can develop. Then 

190 Ufc may lead to a pairing of c-band electrons and /-band 

191 holes, i.e., to the formation of excitons and, provided a 

192 large enough number of excitons is created, a subsequent 

193 spontaneous condensation of these composite quasipar- 

194 tides may develop. In real systems this excitonic insta- 

195 bility is expected to occur, when semimetals with very 

196 small band overlap or semiconductors with very small 

197 band gap are cooled to extremely low temperatures. ^^''^'' 

198 The excitonic condensate then typifies a macroscopic 

199 phase- coherent insulating state, the EI, which separates 

200 the semimetal from the insulator. "^^'"^^ From a theoretical 

201 point of view, Falicov-Kimball-typc models seem to be 

202 the most promising candidates for realizing collective ex- 

203 citon phases. This holds particularly for the generic two- 

204 band extended Falicov-Kimball model (EFKM), where 

205 an EI ground state has been proven to exist in ID and 2D 

206 by constrained-path Monte Carlo simulations. ^^'^'^ Sub- 

207 sequent Hartree-Fock, slave-boson and projector-based 

208 renormalization techniques yield the 2D EFKM ground- 

209 state phase diagram in even quantitative accordance with 

210 unbiased Monte Carlo data,"^^"^^ supporting the applica- 

211 bility of these analytical approaches also in 3D and for 

212 more complicated situations. 



The electronic part of our Hamiltonian, 



H, 



mEFKM 



— He + He 



(7) 



170 Here N denotes the total number of lattice sites. 

171 Regarding the isotropy (anisotropy) of the valence 

172 (conduction) band(s) the Fermi surface of lT-TiSe2 is 

173 only poorly nested, "'^^ which rules out a nesting mech- 

174 anism for the CDW formation even in a simplified 2D 

175 setting. 



214 can be viewed as a multiband extended Falicov-Kimball 

215 model (mEFKM). We note that the mEFKM was studied 

216 previously and has been shown to reproduce the angle- 

217 resolved photoemission spectroscopy data for lT-TiSe2 

218 at temperatures below the critical temperature^'^^'"'^"'*^, 

219 as well as above but close to the critical temperature. ^'^'^^ 

220 We note that the mEFKM exhibits a particular U(l) 

221 symmetry. This can be seen by applying the unitary 



176 2. Electron- electron interaction 

177 Due to the strong screening of the Coulomb interac- 

178 tion in lT-TiSe2,^^ we assume a local electron-electron 

179 interaction. 



222 transformation Uir, 



Ucc 

N 



Cl, I _^C. C|^,aCi^, I 



/ . 7 . / >k+aa'-ka'-k'3^k'-l-a3 
k,k',q a l3>ot 

~1^ 2^ Z^ •/k+q./k'^k'Q'^k'+qQ ' (") 

k,k',q a 






f(t) 



with5„ = iE.(///. 



223 cj^c,-^). The operators fl''' and c\^^ annihilate (create) 



„(t) 



224 an electron at Wannier site i. Obviously we have 



-f^niEFKM — C/^,Q -ff mEFKM C/Ji,Q ■ 



(8) 



225 This symmetry leads to a degeneracy between chiral and 

226 nonchiral CDWs (see below). 

To proceed, we perform a Hartrce-Fock decoupling of 



the electron-electron interaction terms: 



]\[ Z^ A^ A^ '^k+qo'^ka'^k'/gCk'+q/g ^ 
k,k',q a p>a 

^^= mm 4o^kc.'^/3 - NUcc Y^ Yi """/^ 



k a fl^a 



a I3>a 



(9) 



J^ A^ Z_^ /k+q/k'^k'a'^k' 



,'+qa 



k,k',q OL 



Ufc X] "" 51 .^k-^k + UfcTlf J2 cLckc 
a k k,Q 

- ^C^/c'V Zl ''" " Zl ^Q" Z! 4+Q„a/k 
a a k 

-E^QaE/kCk+Q.. + ^|AQap. (10) 



Here we introduced the EI order parameter functions 



A 



■Qa 



Ufc 

N 
Uf 



E(/. 



k'^k+Q„a/ 



^*Qa - "ITT 7 .('^k+Q„a/k)- 



(11) 

(12) 



228 Since the experiments on lT-TiSe2 suggest that the spon- 

229 taneous hybridization of the valence band with one of 

230 the three conduction bands is the dominant effect of the 

231 electron-electron interaction,^ in deriving Eq. (10), we 

232 neglected all terms that mix different conduction bands. 

233 The resulting decoupled Hamiltonian takes the form 

i?mEFKM = E ^k//k/k + E ^kaCko^C^^ 

k k,Q: 

~ Z^ ^Qa'^k+Q^a/k ~ Z^ ^Qa/k'^k+Q«a 
k,a k.Q 

-NUfcUf Y "a - NUcc y^ y^ ngnp 

(13) 



a I3>a 



-j^i:iA«»p: 



235 
236 

237 
238 



with shifted / and c band dispersions: 

Ek/ = ek/ + t//cE'^«' (14) 

a 
£ka = Eka + UfcUf + Ucc E "/^ ' ^^^^ 

The EI low-temperature phase is characterized by non- 
vanishing expectation values (/kCk+Q^„), (Ck+Q„a/k>' 
which cause a correlation gap in the excitation spectrum. 
The mean local electron density in the EI phase is 



Ui 



= 1 + ^ E I (4+Q„./k) I cos(Q„R, + 9c), (16) 

k,a 



where 
1 



k 



QJk 



4ZK4 



k+Q, 



,a/k 



A 



Qa 



Ufc 



240 
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(17) 

Comparing Eq. (16) with relation (1) we recognize the 
amplitude of the charge density modulation as the mod- 
ulus of the hybridization function ^j.(c],,q „/fc)- Like- 
wise we can identify the initial phases 9a in the density 
modulation as the phases of the hybridization functions 
(which coincide with the phases of the EI order parame- 
ters) . 

Note that previous theoretical studies of the 
mEFKM^"'^^ did not include the phase differences of the 
9a, which will be essential for the establishment of a 
chiral CDW.^"® If one is not concerned with the chiral 
CDW problem, disregarding the phases 9a seems to be 
justified since the U(l) symmetry of the mEFKM pre- 
vents the appearance of a stable chiral CDW anyway. 
We show this by analyzing the behavior of the elec- 
tron operators under the unitary transformation U^p^a'- 

hybridization functions (in real space) then transform 

the 



"'^(4/.)e 



iQoRi 



That 



as (cL/,)e'Q"^' = 
phases 9a can be controlled by the unitary transfor- 
mation through the angles ^p. However, in view of (8) 
the total energy is independent of the 9a- Hence these 
phases can be chosen arbitrarily and there is no mecha- 
nism that stabilizes a given phase difference. Therefore 
the mEFKM is insufficient to describe a chiral CDW in 
lT-TiSe2. In the following we will demonstrate that the 
coupling of the electrons to the lattice degrees of free- 
dom can break the U(l)-symmetry of the mEFKM and 
consequently can stabilize a chiral CDW. 



269 B. Lattice degrees of freedom 

270 1. Electron-phonon coupling 

271 For lT-TiSe2 there are experimental and theoretical 

272 evidences that the weak periodic lattice distortion ob- 

273 served comes from a softening of a transverse optical 

274 phonon mode.^®^^*' Of particular importance in this re- 

275 spect is the BZ boundary phonon that connects the F- 

276 and Ljj-points. The valence-band states at the F-point 

277 and the conduction-band states at the La-points are clos- 

278 est to the Fermi energy and, for momentum conserva- 

279 tion reasons, the transfer (excitation) of an electron from 

280 F to La should be accompanied by the creation of a 

281 phonon with momentum — Qq. Hence we assume a sim- 

282 plified single-mode electron-phonon interaction that con- 

283 tains only f-Ca band-mixing terms (cf. Ref. 11): 



H, 



c — ph 



;^EEk(^U^-.)/, 

^ k,q a 

+ (5^")*(&Lq + fcq)4+q./k 



k'^k+qa 



(18) 



284 Here the operator bq describes the annihilation (cre- 

285 ation) of a phonon carrying momentum q, and gf," de- 

286 notes the electron-phonon couphng constant. The Hamil- 

287 tonian (18) can be derived microscopically following the 

288 approach in Ref. 47. Most notably the electron-phonon 

289 coupling (18) breaks the U(l) symmetry of the mEFKM, 

290 i.e., the arbitrariness with respect to the phases 9a is 

291 eliminated. 



2. Phonon-phonon interaction 



316 that the Qi, Q2, and Q3 phonons become soft at the 

317 same temperature;^* we therefore assume I^Q | = {gL \ = 
^'^ Iffo I ~ .9Q (using the notation g^ = .gqe**") and 

319 a;(Qi) = a;(Q2) = ^(Qa) = u. A finite displacement 

320 of the ions is characterized by {bL ) = (&_q ) 7^ 0. We 

321 denote the static lattice distortions by 



Within the harmonic approximation, the Hamiltonian 
for the (nonintcracting) phonons reads"'* 



i^ph = Yl ^^(q)^l^q ' 



(19) 
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296 
297 
298 
299 
300 
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302 



where a;(q) is the bare phonon frequency. A coupling be- 
tween the lattice vibrations results from the anharmonic 
contributions in the expansion of the potential for the 
ions.^^ As we will see below, such an explicit phonon- 
phonon interaction will be necessary to stabilize the chi- 
ral CDW phase. Here we will take into account only 
terms of cubic and quartic order in the lattice displace- 
ment. Then we obtain 



H, 



ph— ph 



' = ^ff Yl -S(qi'q2,q3)(&]ji +&-qj 



N 



q1.q2.q3 

^ibl+b_^J{bl+b_^J 

£'(qi,q2,q3,q4)(6l, 



1 



E 

qi.q2,q3,q4 



qi 



-qi> 
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304 
305 
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309 3. Frozen-phonon approach 

310 We now apply the frozen-phonon approximation and 

311 calculate the lattice distortion at low temperatures. As 

312 elaborated in Refs. 17-20 the phonons causing the lat- 

313 tice displacements in lT-TiSe2 have the momenta Qq, 

314 shown in Fig. 2. Their softening is inherently connected 

315 to strong electronic correlations.^'^ It has been suggested 



5Qa 



Jqq 



2 

~N 
2 

~N 



■9q(^Qo 



-i<E>Q 



^i'^Q 






(25) 

(26) 



322 Obviously, their phases are composed of two contribu- 

323 tions, (j)a — $Q -l- 4>ba, where $„ stems from the phase 

324 of the corresponding electron-phonon coupling constant 

325 and (j)ba originates from the complex average value (6q„ ) . 

326 It will become apparent below that (pta reflects the in- 

327 teraction between the different CDW components. 

If we only take into account the Q^-phonon modes 
the cubic terms in 7fph-ph vanish, since the finite 
^-component of the ordering vectors prevents con- 
straint (23) from being satisfied. Hence the first anhar- 
monic corrections are of quartic order. Here the geometry 
in the xy-pleme is crucial. Replacing all phonon opera- 
tors by their averages, the Hamiltonian H — Hc + Hc-c + 
Hc-ph + ^ph + ffph-ph becomes an effective electronic 
model: 



H =H, 



mEFKM 



DN 



J2 (^Q"4+Q„a/k + ^Qafl 



k^k+Q„Q 



X (^L + ^-q2 ) (bl + 6-q3 ) K, + 6_q J . (20) 

The explicit expressions of i3(qi,q2,q3) and 
^(qi7q2,q3,q4) are lengthy. We only note the 
symmetry relations 

B(-qi,-q2,-q3) = S*(qi,q2,q3), (21) 

-D(-qi,-q2,-q3,-q4) = £'*(qi, q2, q3, q4) , (22) 

and point out the constraints 

S(qi,q2,q3) oc (5qj+q2+q3,G, (23) 

-D(qi,q2,q3,q4) ex (5q,+q2+q3+q,,G ■ (24) 

Here G is a reciprocal lattice vector of the undistorted 
lattice. 



+ (^Q«^Q/3)^e^''(*"+*«) + (<5Q.<5Q,)^e-^''(*-*^) 
+ ('5q.<52j,)V^(*-*^)1+4EI^Q"I' 



TV^ftw 



\Sqo.\' 
4(5q)^ 



(27) 



328 Here D = i?(Qa, Qa, Q/3, Q/3) (/? ^ a), D = 

329 £'(Qq, Qa, ^Qq, ~Qa), and we assume for simplicity 

330 that the phonon-phonon interaction constants are the 

331 same for all combinations of a and 13. D and D are 

332 real. In Eq. (27), the term proportional to D, coming 

333 from the expansion of the phonon-phonon interaction, 

334 guarantees that the free energy is bounded from below 

335 within our approximations. Note that the leading terms 

336 of the phonon-phonon interaction expansion relate the 

337 phases 0^ to each other. Obviously a finite lattice dis- 

338 tortion causes a hybridization between the valence and 

339 the conduction bands. As a consequence a gap in the 

340 electronic spectrum opens, just as in the course of exci- 

341 ton condensation (cf. Eq. (13)). The corresponding local 

342 electron density is given by Eq. (16). 

343 Of particular interest are the phases 9a- Owing to 

344 the second term on the rhs. of Eq. (27), these phases are 



345 coupled to the phases of Squ ■ Let us analyze the possible 

346 values of the phases of the static lattice distortion. We 

347 first note that every Qq is half a reciprocal lattice vector 

348 in the normal phase, i.e., e^'*^°^' = 1, where Ri is a 

349 lattice vector of the undistorted lattice. Hence 



ground-state energy per site is 

N=N^ l'5Q"ll(4+Q„a/k>l COS( 



N N 

4:D 



' OL YOi ) 



k.a 



^Qo 






-iQQRi+2iQoRi 



^-Qo 



(28) 



4 ■ E E l'5Qan'5Q/JpCOs(20„ - 2$„) 
^Q a /3>Q 

cos(20^~2<i>^) + -5-El'^Q"l' 
5q „ 



E^ 



\5qo 



E, 
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That is 5q^ and 6_q create the same phonon. This 

implies (&q^) = (^-q^) = (^Qq)' ^^'^ (^Qo) becomes 
a real number. The phase of SQa is then exclusively 
determined by the corresponding phase of the electron- 
phonon coupling constant $„, cf. Eqs. (25) and (26). 
However, since a triple CDW is not a simple superposi- 
tion of three single CDWs the situation is more subtle. 
Here the change of the periodicity of the lattice caused by 
one CDW component affects the formation of the other 
two components. To elucidate this in some more detail 
let us assume that phonon 1 softens at Tc, while phonon 
2 and phonon 3 soften at Tc — ST. As a result of the 
transition 1 at Tc the periodicity of the crystal changes 
and consequently the BZ changes too (cf. Fig. 4). The 



4(5q)^ 



n iEFKM 

N 



(29) 



376 The first term on the rhs. of Eq. (29) is minimized with 

6l„-0„=7r +s-27r, s = 0,1,2,... (30) 

377 We find that the phases 9a are locked to (j)a + (2s -I- l)7r. 

378 Thus, the phase relationship between (j)i, 02, and (jy^ is 

379 crucial. According to Eq. (29) only the term coming from 

380 the explicit phonon-phonon interaction, i.e., the second 

381 term on the rhs., will affect these phases. Neglecting the 

382 phonon-phonon interaction 0i, (f)2, and 03 are indepen- 

383 dent of each other and chiral and standard CDWs are 

384 degenerate. So the particular form of i?ph~ph is essential 

385 in order to obtain a chiral CDW. 



365 
366 




Ql 



C. Self-consistency equations 

Proceeding further, we assume that the electron- 
phonon interaction constants are real, i.e., $i = $2 = 
$3 = 0. Taking into account the symmetry of the con- 
duction bands and the equality of the electron-phonon 
interactions, we have |i5qi| = |<5q2| = I'^Qsl = I'^qI and 
IAqiI = |Aq2| = |Aq3| = |Aq|. Therewith the free 



393 energy follows as 



FIG. 4. (color online) BZ in the normal phase (black dotted 
hexagon) and (artificial) BZ that would emerge if only the 
phonon Qi softens (filled gray hexagon). Red, green, and 
blue arrows indicate the ordering vectors Qi, Q2, and Q3, 
respectively. 
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vectors Q2 and Q3 are no longer half-reciprocal lattice 

367 vectors and Eq. (28) does not apply. Hence, at Tc ~ ST, 

368 (&Q2) and (fogs) are complex numbers with phases that 

369 have to be determined by minimizing the free energy. For 

370 ir-TiSe2, 6T = 0, but nevertheless the above discussion 

371 remains valid. That is the freedom to fix the phases of 

372 the lattice distortions in an appropriate way results from 

373 the fact that one triple- CDW component must develop 

374 in a lattice structure which is already distorted by the 

375 other two triplc-CDW components. 

Based on the model (27), at zero temperature, we 
can relate the phases of (5qq and J^ki'^k+o^afu) ■ ^^^ 



-^-\Sq\'^ cos(20i) cos(2(/)2) + cos(20i) cos(203) 
5q L 



-i Q n 

+ cos(202)cos(203) +—\SQf, 



(31) 



394 where Uki, = [l + e/^^^"^'^^] . The quasiparticle en- 

395 ergies E^v are obtained by diagonalizing the Hamilton 

396 matrix 



[H] 



f ^kf ^Ql ^Q2 ^Q3 \ 

Aqi Ek+Qii 

Aq2 ek-|-Q22 

VAq3 ek+Q,J 



(32) 



397 Evidently the three conduction bands arc coupled to 

398 each other via their coupling to the single valence band. 

399 The resulting quasiparticle band structure exhibits a gap, 

400 whose magnitude is determined by the squared modulus 

401 of the "gap parameter" 



A 



Qq 



5Qq 



^Qq- 
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Naturally the relation (30) maximizes the gap parameter 
|Aqq| for every given |(5qq,| and |Aqq|. On the other 
hand, i?kiy depends, via cos{9a — 4'a)j on the phase dif- 
ferences only, i.e., the quasiparticle band structure can 
not be used to fix the particular values of 9i, 02, and 63. 
We further note that |Aq| > if and only if \Sq\ > 0. To 
elaborate let us first consider the case gq = 0. In the EI 
phase (|Aq| > 0) the system realizes a CDW. When ^q 
becomes finite in addition, the lattice adjusts commen- 
surate with the electron density modulation. Hence, in 
this case, any finite ^q immediately results in |5q| > 0. 
On the contrary, at vanishing Coulomb interaction but 
sufficiently large gq > 5q,c, a lattice instability devel- 
ops leading to a finite X^ki^k-t-QaQ'^k)- "^^^^ hybridiza- 
tion parameter enters the explicit equation for the EI 
order parameter (12). |Aq| > then follows from any 
finite Coulomb interaction. Our approach therefore does 
not discriminate between an excitonic and phonon-driven 
instability if both electron-electron and electron-phonon 
interactions are at play. 

We now relate |Aq| and \Sq\. Minimizing F, Eq. (31), 
with respect to 9a yields the constraint (30). This entails 
that the gap parameter is given by |Aqq| = Ii^qI + |Aq|. 
Depending on the sign of D, a different phase relationship 
between (j>i, <f>2, and 03 will be favored. For D > we 
find that one phase has to be a multiple of n and another 
phase must be an odd multiple of tt/2. The third phase 
is arbitrary. In this case the considered phonon-phonon 
interaction introduces a fixed phase difference between 
the lattice distortion modes. That is, 9i 7^ 6*2 7^ ^3, and 
we observe a chiral CDW. On the other hand, for D < 0, 
F takes its minimum if all phases are a multiple tt/2. The 
resulting CDW is then nonchiral. For both cases (D > 
and D < 0), the free energy is minimized at finite values 
of \Sq\ if the constraint 



D> A\D\ 



Q 9q 

+ COs(2(^i) COS(203) + C0S(2(^2) cos(2(/.3) 



9q 



Let us finally stress that we can only specify values 
for the phases 9a and 4>a- Which particular phase takes 
one of these values remains open. For instance, the si- 
multaneous transformations 02 -^ 9^ and 02 -^ 03 do 
not change the free energy, but convert a clockwise chi- 
ral CDW in an anticlockwise one. The degeneracy of 
these two CDW states is in accord with the experimental 
findings for lT-TiSe2.^ 



D. Characterization of the CDW state in ir-TiSe2 

Experiments identify a close connection between the 
appearance of the CDW state and the periodic lattice dis- 
placement in lT-TiSe2.^ In the scheme developed in the 
preceding sections the CDW is directly accompanied by 
a lattice distortion. There the j-component {j = x, y, z) 
of the lattice deformation caused by the phonon a is 



<(r) 



1 



E 



2h 



K^Q )|e;,j(Qa) cos(Qar-0(,a), 



(36) 

where e/j(QQ,) is the j-component of the polarization 
vector, I labels the atoms in the unit cell, and Mi is the 
mass of the atom I. Clearly each CDW component ol 
produces a 3D lattice distortion. Since 061 7^ 062 7^ 063 
if 01 7^ 02 7^ 03 , the lattice will be differently affected by 
the phonons Qi, Q2, and Q3. Of course the lattice defor- 
mation by the phonon mode Qq, is position-dependent; in 
this way a complicated distortion pattern of the ions can 
occur. An instructive picture can be achieved, however, 
if one neglects the r-dependence in the a;j/-plane. In this 
simplified situation, depending on the z-componcnt as a 
function of r, the magnitude of the lattice displacement 
differs along Qi, Q2, and Q3. As a result the differ- 
ent ionic layers of lT-TiSe2 are dominated by different 
phonon modes. ^ The situation where the lower Se-ion 
layer is largely affected by the phonon mode Q3, the Ti- 
ion layer by phonons with momentum Q2, and the up- 
per Se-ion layer by the Qi phonon mode, is illustrated 
schematically in Fig. 5(a). Let us consider the upper 



(34) 



437 is fulfilled. Since the experiments suggest a chiral CDW, 

438 we take Z) > 0. Without loss of generality we further set 

439 01 = 0, 02 = 7r/2, 03 = 37r/10, and assume D = 413. 

440 The moduli! of |(5q| and |Aq| arc then determined by 

441 dFld\bci\ = and 5F/9|Aq| = 0, respectively. Straight- 

442 forward calculation gives 

-^|Aq| ^ ^I^qI +84l'5Qpfcos(20i)cos(202) 

Ufc 45, 




FIG. 5. (color online) (a) For a chiral ordering the maximum 
lattice distortion due to the phonon Qq may be located in 
different ionic layers, (b) The Se ions in the upper layer are 
differently affected by the phonons having momentum Qi, 
Q2, or Q3. For further discussion see text. 
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479 plane of Sc ions, which is analyzed in scanning tunneling 

480 microscopy experiments.^ There a relative difference of 



481 the phases (j)ba leads, e.g., to a stronger displacement of 

482 the ions in the direction of Qi than in the direction of Q2 

483 and Q3 [see Fig. 5(b)]. Then the CDW transition can be 

484 viewed as a formation of "virtual layers" with ordering 

485 vectors assigned to a helical structure.^ The described 

486 distortion scenario equates with a fixed phase difference. 

487 Thereby the only crucial parameters arc 0i, 02, and ^3; 

488 the finite z-componcnt of the ordering vectors is not a 

489 required prerequisite for the chiral CDW. 

490 We add that in our model both CDW transitions ap- 

491 pear simultaneously, i.e., Tcdw = T'chirai- This contra- 

492 diets the result Tchirai < ^cdw by van Wezel.^'^ In that 

493 Landau-Ginzburg approach terms of second and quartic 

494 order in the order parameter appear. While the second- 

495 order term favors equal phases, the quartic-order term 

496 supports a phase difference.^ Hence there is a compcti- 

497 tion between a nonchiral and a chiral CDW state. This 

498 competition is absent in our approach where we have 

499 only a single term that fixes the values of 0i, 4>2, i's re- 

500 spectively ^i, 02, Oz- Very recent X-ray diffraction and 

501 electrical transport measurements provide some evidence 

502 that the transition to the chiral CDW takes place at lower 

503 temperature than the CDW transition. ^° The difference 

504 between Tcdw and Tchirai is less than ten Kelvin how- 

505 ever. Having disregarded fluctuation corrections, which 

506 are known to play an important role in the vicinity of 

507 a phase transition, we can not conclusively resolve this 

508 issue. 



III. NUMERICAL RESULTS 



A. Model assumptions 

In view of the quasi-2D crystallographic and electronic 
structure of lT-TiSe2, and in order to simplify the nu- 
merics, we restrict the following analysis to a strictly 2D 
setting. Moreover, being close to the Fermi energy, we 
will approximate the bands parabolically^° 



£k/ 



-tf ki 



rkl = t'cikx - Qlxf + tliky - Qly)2 + E, 



(37) 
(38) 



516 with hopping amplitudes tf, t^, and t^. The other two 

517 conduction bands ek2 and Sks have analogous dispersions, 

518 but the momenta are rotated by 27r/3 and 47r/3, respec- 

519 tively. All three conduction bands share the same mini- 

520 mum Ec, see Fig. 3. 

521 From the band dispersion provided by Monney et at in 

522 Rcf. 10 we derive tf = 1.3 cV, which will be taken as the 

523 unit of energy hereafter, and t^ = 0.042 and t^ = 0.105. 

524 The bare phonon frequency is estimated as hu ~ 0.013, 

525 in accordance with the value given by Weber et al. in 

526 Ref. 20. Furthermore, we set Ec = —3.30 and Ucc = 

527 Ufc + 1.0. For the phonon- phonon interaction constant 

528 we take D = 10~^. 

529 The self-consistency loop, comprising the detcrmina- 

530 tion of the EI order parameter, the static lattice dis- 

531 tortion, the total and partial particle densities and the 



chemical potential, is solved iteratively until the relative 
error of each physical quantity is less than 10~^. The nu- 
merical integrations were performed using the Cubpack 
package. ^^ 



536 B. Phase diagram of the mEFKM 

537 To set the stage for the analysis of the interplay of 

538 Coulomb and electron-phonon interaction effects we first 

539 discuss the phase diagram of the pure mEFKM, cf. Fig. 6. 

540 Here, since gq = (and as a result (5qq = 0), the EI 

541 low-temperature phase typifies a normal CDW. As for 

542 the EFKM on a square lattice (see inset), at T = we 

543 find a finite critical Coulomb strength above which the 

544 EI phase does not exist. This is because the large band 

545 splitting caused by the Hartree term of the Coulomb in- 

546 teraction prevents c-f electron coherence. ^^ In contrast 




FIG. 6. Phase diagram of the mEFKM. The inset displays the 
schematic phase diagram of the simplified two-band EFKM 
on a square lattice according to Ref. 42. 



548 
549 



to the EFKM, in our four-band model we also find a 

550 critical lower Coulomb strength for the EI phase. This 

551 can be understood as follows: since the valence band 

552 is isotropic while the conduction-band dispersions are 

553 strongly anisotropic, particles close to the Fermi level 

554 do not find a large number of partners with appropri- 

555 ate momentum for electron-hole pairing. Thus, for Ufc 

556 smaller than a critical Coulomb attraction, the amount 

557 of energy to create a macroscopic number of excitons is 

558 larger than the energy gain from the condensation tran- 

559 sition into the EI state. Therefore the system remains 

560 in the semimetallic phase. ^^ The rather abrupt increase 

561 of the critical temperature at the lower critical Coulomb 

562 interaction is due to the degeneracy of of the conduction 

563 bands and the particular anisotropy used. 



564 C. Influence of the lattice degrees of freedom 

565 We now analyze the situation when phonons partic- 

566 ipatc in the CDW formation. In Fig. 7 we consider 

567 the dependence of the EI order parameter |Aq| and the 

568 static lattice distortion \5q\ on the interaction constants 
5B9 at T = 0. For small electron-phonon coupling the be- 

571 havior of the EI order parameter resembles the situa- 

572 tion for the mEFKM (EI limit). Differences appear for 
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FIG. 7. EI order parameter |Aq| (left panel) and static lattice 
distortion |5q| (right panel) as a function of the Coulomb 
interaction C//c at T = 0. Results are given for gq = 0.02 
(solid lines) and gq = 0.11 (dashed lines). For comparison, 
in the left-hand panel, the EI order parameter of the pure 
mEFKM is included ((7q = 0; dotted line). 



573 gq = 0.11 where we do not find a lower critical C/yc for 

574 the CDW phase, i.e., now any finite Ufa leads to a finite 

575 |Aq|. The order parameter |Aq| shows an almost linear 

576 dependence on V fci which means that the electron- hole 

577 pairing probability (oc |Aq|/[//c) remains approximately 

578 constant. Tlic tendency towards cxciton formation and 

579 condensation is notably enhanced as gq increases at fixed 

580 V fc- The static lattice distortion is tiny and follows the 

581 behavior of the EI order parameter for small gq. For 

582 (7q = 0.11 we observe a non-monotonic behavior of |(5q| 

583 as a function of the electron-electron interaction. At 

584 very large electron-phonon couplings the infiuence of the 

585 electron-electron interaction on |(5q| weakens (not shown 

586 in Fig. 7). 

587 Figure 8 gives the dependence of the order parame- 

588 tors on the electron-phonon coupling constant at fixed 

589 C//c. For V fc = 2.5 the order parameters remain finite 

591 as gg — > since we already have a CDW in the pure 

592 mEFKM (sec Fig. 6). A Coulomb strength C//c = 5.0, on 

593 the other hand, does not lead to an EI if phonons are ab- 

594 sent. Then wc find a lower critical gq. Furthermore we 

595 observe a saturation of the EI order parameter with in- 

596 creasing gq, since |Aq| is limited by the number of elec- 

597 trons that can participate in the electron-hole pairing. 

598 Within our approximation the expectation value of the 

599 phonon-creation operator is not bounded and |(5q| con- 

600 tinuously increases with enhanced electron-phonon cou- 

601 pling. Note that an increasing phonon-phonon inter- 

602 action strength D reduces the static lattice distortion 

603 and therefore leads to a smaller EI order parameter (not 

604 shown). 

605 Finite temperature effects are illustrated in Fig. 9. For 

606 very small electron-phonon coupling the phase diagram 
608 again resembles the situation for the mEFKM. As the 
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FIG. 8. EI order parameter |Aqt (left panel) and static lattice 
distortion |(5q| (right panel) as a function of the electron- 
phonon coupling (/q at r = 0. Results are given for [7/^ = 2.5 
(solid lines) and f//c = 5.0 (dashed lines). The right-hand 
inset enlarges the region of small gq. 




FIG. 9. CDW phase boundaries for yq = 0.02 (solid hue) 
and gq = 0.11 (dashed line). The inset gives the EI order 
parameter (solid line) and the static lattice distortion (dashed 
line) as a function of temperature for gq = 0.11 at U jc = 2.5. 



609 interaction strength yq increases the situation changes 

610 dramatically. For sufficiently large electron-phonon cou- 

611 plings, we no longer find critical lower and upper val- 

612 ues U fc for the CDW transition and the transition tem- 

613 perature increases linearly with V jc- That is the criti- 

614 cal temperature is significantly enhanced by (7q. Obvi- 

615 ously electron-hole attraction and electron-phonon cou- 

616 pling act together in the formation of a very stable CDW 

617 phase. The inset shows how the order parameters are 

618 suppressed by increasing temperature. Here the static 

619 lattice distortion reveals a typical mean-field dependence. 

620 Our findings arc summarized by the ground-state 

621 phase diagram shown in Fig. 10. For weak clectron- 

622 phonon couplings (/q the CDW is mainly driven by the 



10 



623 Coulomb attraction Ufc between electrons and holes. 

624 The greater gq, the larger the region where the CDW 

625 is stable. For gq > 0.09 the electron-phonon coupling 

626 alone can cause the CDW transition, even at Ufc = 

627 (blue line in Fig. 10). The CDW is chiral in this limit, 

628 whereas the CDW in the opposite EI limit does not ex- 

629 hibit chirality (gq = 0, red line in Fig. 10). 
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FIG. 10. (color online) Ground-state phase diagram of the 
niEFKM with additional electron-phonon coupling. The 
CDW phase is characterized by a finite gap parameter |Aq|. 
The red line at gq = marks the EI phase of the pure 
mEFKM. The blue line at Ufc = refers to a CDW in- 
duced solely by the electron-lattice interaction. The green 
point designates the range of model parameters appropriate 
for ir-TiSe2. 



632 D. Relation to lT-TiSe2 

633 Based on the phase diagram derived for the mEFKM 

634 with electron-phonon coupling, we now attempt to cs- 

635 timate the electron-electron and electron-phonon intcr- 

636 action constants, Ufc and gq, for lT-TiSe2. To make 

637 contact with the experiments we take the displacements 

638 of the Ti ions measured by Di Salvo et al.: ui=ti = 0.04A 

639 with u = 'Ylii ^i-'^ Then, from Eq. (36), we can specify the 

640 value of |(foq)|. For lT-TiSc2, the gap parameter was de- 

641 termined experimentally as 120 meV by Monney et al., 

642 cf. Ref. 44. Adjusting this value to our theoretical re- 

643 suits yields Ufc « 2.5 (« 3 eV) and 5q « 0.03 (« 0.04 

644 eV); see the green marker in Fig. 10. For these values 

645 both the theoretical ion displacement and gap parame- 

646 ter are in the same order of magnitude as the measured 

647 ones. Using Ufc ~ 2.5 for lT-TiSe2, the electron-hole 

648 pairing is BCS-likc. Since yq ~ 0.03 is too small to 

649 cause a CDW for vanishing Coulomb interaction and, as 

650 discussed above, the EI scenario alone will not yield a 

651 stable chiral CDW, our results are in favor of a combined 

652 lattice-deformation/EI mechanism for the experimentally 

653 observed chiral CDW transition, as suggested in Rcfs. 22 



and 53. 



655 IV. CONCLUSIONS 

656 In this work we have argued how the chiral charge- 

657 density-wave (CDW) phase in lT-TiSe2 may be sta- 

658 bilized. In the framework of the multiband extended 

659 Falicov-Kimball model (mEFKM) we showed that a 

660 purely electronic — exciton pairing and condensation — 

661 mechanism is insufficient to induce the observed (long- 

662 ranged) chiral charge order. We propose that the cou- 

663 pling of the electrons to the lattice degrees of freedom is 

664 essential for the formation of a chiral CDW state. In par- 

665 ticular, an explicit phonon-phonon interaction between 

666 the phonons that become soft at the critical temperature 

667 has to be taken into account in order to stabilize the chi- 

668 ral CDW. We note that in our model clockwise and an- 

669 ticlockwisc CDWs arc degenerate. This is in accord with 

670 experimental findings.^ The chiral property can properly 

671 be observed in the ionic displacements accompanying the 

672 CDW in ir-TiSe2. 

673 Concerning the microscopic mechanism underlying the 

674 CDW transition, we demonstrated that electron-electron 

675 interaction and electron-phonon coupling support each 

676 other in driving the electron-hole pairing and finally the 

677 instability. This suggests that the CDW transition in 

678 lT-TiSe2 is due to a combined lattice distortion and 

679 exciton-condensation effect. The outcome is a spon- 

680 taneous broken-symmetry CDW low-temperature state 

681 with small but finite lattice deformation. Of course, both 

682 the mean-field treatment of the Coulomb interaction and 

683 the frozcn-phonon approach are rather crude approxima- 

684 tions and a more elaborated study of the complex in- 

685 terplay between the electronic and phononic degrees of 

686 freedom is highly desirable to confirnr our proposed scc- 

687 nario for the chiral CDW transition in ir-TiSe2. 

688 Let us finally point out that we called |Aq| the 

689 excitonic- insulator order parameter on account of its ana- 

690 log in the EFKM.35-42 ^^j^g meaning of a finite |Aq| in 

691 the presence of a band coupling is imprecise however. 

692 Likewise a spontaneous hybridization of the valence band 

693 with one of the conduction bands, signaling the exciton 

694 condensate in the mEFKM, may be induced by a suf- 

695 ficiently large electron-phonon coupling. A general cri- 

696 terion for the formation of an exciton condensate in a 

697 strongly coupled band situation has not been established 

698 so far. This is an open issue which deserves further anal- 

699 ysis because of its relevance in characterizing the nature 

700 of CDW transitions also in other materials. ^''^'^^ 
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